A deterministic origin for instabilities in a magneto-optical trap 
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The paper reports on the observation of a new type of instabilities in the atomic cloud of a 
magneto-optical trap (MOT). They appear as large amplitude self-oscillations, either periodic or 
erratic, depending on the parameters. An original model taking into account the shadow effect and 
an homogeneous distribution of the atoms in the cloud is developped. It allows us to show that 
these instabilities originate on the one hand from an anomalous Hopf bifurcation, and on the other 
hand from the extreme sensitivity of the system to noise. 
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The magneto-optical cooling of atoms is at the origin of 
a renew of the atomic physics. It is used in various fields, 
as the realization of Bose-Einstein condensates |[ or the 
study of diffusion in lattices Q , and could lead to several 
applications, as atomic clocks or quantum computing 
H . Although the technology and realization of magneto- 
optical traps (MOT) is well mastered, some experimental 
adjustments remain empirical. It is in particular well 
known that for dense atomic clouds close from resonance, 
instabilities appear in the spatio-temporal distribution of 
the atoms. This problem is usually fixed by empirically 
misaligning the trapping beams. 

A recent study has shown that instabilities can have a 
stochastic origin, through a phenomenon of coherent res- 
onance amplifying the technical noise [0 . It also showed 
that the shadow effect is at the origin of these instabili- 
ties : because of the absorption of light inside the cloud, 
the intensities of the backward and forward beams are lo- 
cally different, leading to an internal attracting force. In 
the configuration where each backward beam is obtained 
by retro-reflection of the forward beam, backward and 
forward intensities are globally different, and an external 
force appears, displacing the cloud along the bisectors of 
the trap beams. 

In the present work, we show that the atomic cloud 
in the MOT can also exhibit deterministic periodic and 
chaotic- like instabilities. The basic model developed in 
f| is not able to reproduce them, and we propose a mod- 
ified model where the cloud is no more considered as 
pointlike. This new model predicts the same type of 
stationary solutions as in S, leading to the same type 
of stochastic instabilities. Moreover, it allows us to re- 
trieve the deterministic instabilities: in particular, it ex- 
hibits large amplitude periodic oscillations induced by 
an anomalous Hopf bifurcation, and a large amplitude 
stochastic behavior. It also establish definitively the role 
of the shadow effect in the dynamics of the atomic cloud 
in the MOT. 

The initial hypotheses are the same as in ||, i.e. a 



ID model taking into account the shadow effect. The 
system is modelized on the one hand through the motion 
equations of the center of mass z of the cloud, and on 
the other hand through a rate equation of the number of 
atoms n in the cloud. We have: 
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where M is the mass of the cloud, Ft the total exter- 
nal force, n e the atom number at equilibrium and B 
the population relaxation rate. As in [5|, n e is assumed 
to depend on z, to take into account the depopulation 
of the cloud when it moves away from the trap center. 
We define a distance zq beyond which the trap is empty 
(n e = 0). For z < zq, we assume a quadratic behavior 

n e = riQ ^1 — (z/zo) 2 ^ where no is the cloud population 

at the trap center ||. 

The innovation of this model concerns the way the 
shadow effect is treated. We consider that, starting from 
an input forward intensity Zi, the intensity after a first 
crossing of the cloud (i.e. the input backward intensity) is 

12 < Ii , and the remaining intensity after a second cross- 
ing of the cloud (i.e. the output backward intensity) is 

13 < I 2 - The rate of photons absorbed in the forward 
(resp. backward) beam is A (Ii — I 2 ) (resp. A (I 2 — h)), 
with A = S/hv, where S is the section of the cloud and hv 
the energy of a photon. The forces associated with each 
beam is the product of the number of absorbed photons 
by the elementary momentum hk. We obtain finally: 
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To get a relation between I±, I 2l and 7 3 , we need 
to solve the equations of propagation of the two beams 
through the atomic cloud. For the sake of simplicity, we 
assume that the atomic transition isaJ = 0^J=l 
transition. As the MOT is operated with high intensity 
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beams and small detunings, a Doppler model is suited. 
Inside the cloud, the intensity I + (7_) of the a + (cr_) po- 
larized forward (backward) beam evolves due to photon 
scattering. The scattering rate for a± photons is propor- 
tional to the corresponding excited state populations H± . 
The evolution equations of the intensity simply writes 
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where T is the natural width of the atomic transition and 
p is the atomic density in the cloud. The populations Tl± 
are given by the steady state of the master equation. The 
underlying hypothesis is that the evolution of the exter- 
nal degrees of freedom is much slower than that of the 
internal ones. The populations II± depend both on I + 
and J_, so that (|J) is a set of coupled nonlinear equa- 
tions. They are integrated numerically from the side of 
the cloud where I + = J_ = h, to the other side, where 
7_ = I3 and I + = Ji, assuming that the density p is 
constant, as a result of multiple scattering |J. 

The control parameters of this system are the same as 
in [||. The most suitable one is the detuning Ao, ex- 
pressed in units of F, between the trap laser beams and 
the atomic transition. The global organization of sta- 
tionary solutions in the parameter space is also the same 
as in [^j. In particular, stationary solutions as a func- 
tion of Ao evolve from a flat dependence to bistability 
when another parameter (e.g. no) is changed. The fold 
appearing for the intermediate values is at the origin of 
the stochastic instabilities reported in 0. 

We focus here on the unstable zone appearing between 
the area studied in j^] and bistability. In this situation, 
the stationary solutions are unstable on the fold (fig. [I]). 
For detuning smaller than the fold, at the left of point Hi 
on fig. [l]a, the fixed point is a stable focus (F): the sta- 
tionary solutions are stable and associated with an eigen- 
frequency ujf decreasing with the detuning (fig. 0b). At 
the edge of the fold, the system exhibits a Hopf bifurca- 
tion (point Hi): the fixed point becomes a saddle- focus 
(SF), and the stationary solutions become unstable. As 
the detuning is further increased, the eigenvalues become 
real in point Pi (fig. pL,b), so that up disappears and the 
fixed point becomes a saddle node (SN). Finally, when 
the detuning is still increased, the inverse sequence ap- 
pears for the fixed point (SN — > SF — > Hopf bifurcation 
-F). 

Let us now examine the dynamical behavior of the 
atomic cloud in the different regions. In the F zone, 
the stationary solution is stable. As the detuning is in- 
creased, the stationary solution becomes unstable in Hi, 
and a stable periodic orbit Cjj appears in the vicinity 
of the fixed point, as it is usual with a Hopf bifurcation. 
The atomic cloud exhibits weak amplitude oscillations 
around the unstable fixed point, with a frequency close 
to Wf. The orbit amplitude increases slowly with the de- 
tuning, but remains very weak (less than 10 nm in z for 
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1. Theoretical evolution of the behavior of the cloud 
as a function of the detuning. In (a), the stationary solution 
z s of z is stable (full line) or unstable (dashed line) . In points 
Hi and H2, a Hopf bifurcation occurs, while in points Pi and 
P2, ujf vanishes. F (focus), SF (saddle focus) and SN (saddle 
node) refer to the nature of the fixed point representing the 
stationary solution in the phase space. In (b), evolution of 
ujf versus Aq. In (c), plot of the instability frequencies uja 
(circles) and ujb (squares). Parameters for the calculations 
are h = 33 mW/cm 2 , p = 2 x 10 10 cm" 3 , n = 6 x 10 8 , 
zq — 3 cm, B — 5 s _1 and a Zeeman shift of 3F cm" 1 . 



the parameters of fig. |l|). 

In fact, if Ao is further increased, Ch disappears for a 
value A^i of Ao very close from Hi. It is replaced by an- 
other limit cycle Ca , which differs from C'h on three main 
points: its amplitude is much larger than that of Ch, its 
frequency uj\ is not u>p, and its shape is not sinusoidal 
at all. Fig. ||a shows an example of a Ca cycle: in each 
period, a slow increasing of z is followed by a fast growth 
where z becomes much larger than z s , and a fast de- 
creasing after which the cycle starts again. Meanwhile, n 
increases slowly and regularly, then decreases. Note that 
in the example of fig. ^|a, the stationary value n s of the 
population is outside the limit cycle, and so the interpre- 
tation of the behavior in terms of the relative values of z 
and n with respect to z s and n s is not straightforward. 
In particular, it cannot be interpreted as a relaxation os- 
cillation process around the unstable stationary values. 
In fact, as we show below, this behavior is not connected 
with the unstable fixed point properties. 

Because Ch becomes unstable for a value so close from 
Hi, the Hopf bifurcation appears as anomalous: while 
such a bifurcation leads usually to progressively growing 
cycle, it generates here a large amplitude periodic orbit. 
However, the Ca behavior is not linked to the SF zone: it 
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FIG. 2. Examples of the behavior of the cloud. The full 
(dashed) line curve is a plot of z (n) versus time. The hor- 
izontal full (dashed) line marks the stationary value z 3 (n s ). 
In (a), n a /n = 0.757 is outside the figure, (a) shows a Ca 
instability for Ao = —0.37 ; (b) shows a Cb instability for 
Ao = —0.35; (c) corresponds to the same parameters as (b), 
but a noise level of 7% has been added on 7i . Other parame- 
ters are the same as in fig. nl. 



still exists beyond Pi and P2, without any discontinuity. 
The amplitude of Ca is several millimeters when it ap- 
pears in Aai, and increases regularly with the detuning, 
so that for a value Aa2 of Ao, located between P2 and H2 
for the parameters of fig. [l], the cloud border reaches Zo. 
At this point, the shape of the limit cycle qualitatively 
changes. Indeed, the atoms beyond zq are lost, and so 
the population in the cloud can decrease rapidly. The re- 
sulting temporal behavior is still a periodic cycle Cb and 
may be described as previously, except that the decreas- 
ing of z is much faster and that of the population much 
larger (fig. ||b). Note that this behavior is observed be- 
yond H2: this means that generalized bistability occurs 
between Cb and the stable stationary solution. This con- 
firms that the periodic instabilities are not directly linked 
to the fixed point properties. 

The Ch, Ca and Cb behaviors differ also by their fre- 
quency. As discussed above, the Ch frequency is up, 
and so varies rapidly as a function of the detuning. The 
frequency uja of the Ca regime is very small when the 
regime appears in Aai (fig- 0c)- When the detuning 
is increased, uja reaches rapidly a value of the order of 



2tt x 10 s _1 , and then remains almost constant. Finally, 
the frequency ujb of the Cb regime is also almost con- 
stant when Ao is changed, with a frequency smaller than 
uja (about in fig. |l]). 

Previous studies have shown the fundamental role of 
noise in this system . The influence of noise on deter- 
ministic instabilities is well known: fixed points and limit 
cycles are usually robust with respect to noise, whose 
main effect is to slightly shift the bifurcation points [Q. 
So we do not expect to observe spectacular changes in 
the stationary, Ch and CU behaviors if noise is added in 
the system, and this is confirmed by the simulations. The 
Cb behavior is different, as the cloud could be very sen- 
sitive to noise in the vicinity of zq: indeed, noise should 
induce large variations in the decreasing of n, and hence 
in the period of the dynamics. This is confirmed by the 
numerical simulations: fig. |^c shows the behavior of the 
cloud for the same conditions as in fig. ||b, except that 
noise has been added on the trap intensity. As expected, 
the dynamics is no more periodic, exhibiting large fluc- 
tuations in the return time of the dynamics. 

To check the existence of the deterministic instabilities 
predicted by the model, we used the same experimen- 
tal setup as described in ||], i.e. a three-arm <r + — <r~ 
magneto-optical trap (MOT) on Cesium, with mirrors to 
produce the counter-propagating waves of each arm, and 
a crossed couple of 4-quadrant photodiodes to monitor 
the cloud. As in ||, forward and backward beams are 
carefully aligned. The only changes with respect to the 
experiment described in concern the parameter values, 
and in particular the intensity of the trap laser beams, 
adjusted to obtain deterministic instabilities. 

As expected, the Ch cycle is not observed. The most 
commonly observed periodic instabilities look as C'a in- 
stabilities (fig. H). Fig. H shows the evolution of the sig- 
nal frequency uj as a function of the detuning. As in the 
model, uj is almost constant on the interval where Ca in- 
stabilities appear, i.e. for —1.7 < Ao < —0.8. The value 
of uj depends on the experimental parameters, as the trap 
beam or the repump intensities. It is typically of the or- 
der of one or a few hertz, in good agreement with the 
model. On the contrary, the detuning interval on which 
instabilities appear is one order of magnitude larger than 
that predicted. This is the main discrepancy between 
the experimental observations and the model. However, 
to make a real comparison, we should take into account 
the inevitable variation of no when Ao is changed. Note 
that in the model, a simultaneous change of no and Ao 
leads to a relative stretch of the unstable zone. Unfor- 
tunately, as we have no way to establish experimentally 
the relation between no and Ao , we are not able to check 
the amplitude of the correction in the present model. 

When the detuning is increased, the behavior shifts to- 
wards Cb instabilities. Because of noise, the transition 
is not abrupt as in the simulations. The limit cycle be- 
comes more and more noisy as Aq is increased, while the 
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FIG. 3. Experimental record of a Ca periodic instability. 
Parameters are h = 11 mW/cm 2 and A = -1.4. The beam 
waist of the trap beams is 5 mm, while the magnetic field 
gradient is 13 G/cm. £ is proportionnal to the mean size of 
the cloud. 



15- 
10- 



n 1 1 1 1 1 1 1 

-2.0 -1.6 -1.2 -0.8 

A (r units) 

FIG. 4. Evolution of the instability frequency u) versus the 
detuning. Parameters are the same as in Fig. ^ except I\ = 20 
mW/cm 2 . 



In conclusion, we have demonstrated the existence of a 
new type of instabilities in the MOT clouds. These insta- 
bilities, contrary to those already observed, are determin- 
istic, implying that a simple amelioration of the experi- 
mental noise cannot improve the cloud stability. Instabil- 
ities consist in large amplitude periodic or erratic cycles. 
We have proposed a new model taking into account the 
intensity distribution inside the cloud. This model allows 
us to interprete the instabilities through an anomalous 
Hopf bifurcation generating the large amplitude pulses, 
and shows that even in the case of deterministic instabil- 
ities, noise can play a main role in the behavior. These 
results are obtained for a three-beam configuration, but 
we are expecting them to be more general. In particular, 
in the six beam configuration, the shadow effect also in- 
duces nonlinearities, leading to similar dynamics on the 
atomic density. However, the present model is still too 
simple to give a quantitative agreement with the experi- 
ments. A 3D model with a non constant atomic density 
will give a finer description of the dynamics. 
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frequency decreases. Finally, for A > —0.55, the insta- 
bilities disappear and the behavior is again stationary. 
As expected in the above discussion, the Cb instabilities 
appear in their erratic form (fig. |^). The main char- 
acteristics of the experimental behavior corresponds to 
the theoretical results: the decreasing stages of z and n 
are very fast, while the mean frequency of the signal is 
shorter than that of Ca ■ 
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FIG. 5. Experimental record of a C2-like instability. Pa- 
rameters are the same as in Fig. W with Ao = —0.6. 
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